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Abstract We study temperature chaos in a two-dimensional Ising spin glass with random 
quenched bimodal couplings, by an exact computation of the partition functions on large systems. 
We study two temperature correlators from the total free energy and from the domain wall free en- 
ergy: in the second case we detect a chaotic behavior. We determine and discuss the chaos exponent 
and the fractal dimension of the domain walls. 
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Introduction A characteristic feature of spin 
glasses is the presence of chaos under small changes in the 
quenched couplings, in the temperature, or in the mag- 
netic field HSU SHH, HUB- With th<5 expression 
temperature chaos we refer to the fragility of the equilib- 
rium states of a disordered system under small temper- 
ature changes. Let us consider two typical equilibrium 
configurations of such a system under the same realiza- 
tion of the quenched disorder: the first configuration is in 
equilibrium at temperature T, while the second one is in 
equilibrium at temperature T 1 = T + AT. One says that 
there is temperature chaos if for arbitrarily small (but 
non-zero) values of AT, the typical overlap of two con- 
figurations at T and T' goes to zero when the system size 
diverges. The spatial distance £(T, AT) over which such 
overlaps decay is called the chaos length, and, as we will 
discuss better in the following, it scales as I ~ AT^ 1 ^ 
when AT -> 0. 

The droplet picture 0, of spin glasses predicts that 

C = y-0, (1) 

where d s is the fractal dimension of the droplet inter- 
faces and 9 is the usual spin glass stiffness exponent. In 
two-dimensional (2D) spin glasses with continuous dis- 
tributions of spin-spin couplings, one has d s w 1.27 (H3) 
and 8 ~ —0.285, leading to the prediction £ w 0.92; di- 
rect measurements of £ 0,0 are in good agreement with 
this prediction. The situation in the model with binary 
couplings (Jij — ±1) is notably different: no study of 
temperature chaos has been performed, but if one works 
exactly at T = using the estimates d^ =0 » 1.0 (flflj) 
and 9 T=0 = (El E3), one would expect C T=0 « 0.5 
which is very different from the continuous distribution 
case, suggesting two universality classes. However, even 
though the T = properties of spin glasses with continu- 
ous and discrete distributions arc different, a new picture 
has recently emerged EH m which one shows that only 
one universality class exists for T > critical properties 



when we consider the limit T — > 0. In our work here we 
measure £ T ^° in the ± J model and show that its value 
is compatible with the one of the model with Gaussian 
couplings, giving further credence to the generalized uni- 
versality scenario 0, . 

Our study is based on exact partition function com- 
putations for large but finite size systems, for many re- 
alizations of the quenched disordered couplings 0, 0] . 
Availability of the full density of states of a large system 
is a very powerful tool and it allows us to investigate the 
chaotic behavior of the free energy for the ±J model. We 
also extract entropies of domain walls and thus provide 
an improved estimate of the exponent d s . 

The outline of the paper is as follows. First we discuss 
the model, our observables and the methods we use. We 
discuss how to detect temperature chaos. We analyze 
what happens for the total free energy, finding there is 
no temperature chaos there. Then we move on to domain 
wall free energies. There we compute the two exponents, 
d s from the entropy of T = domain walls, and £ from 
T > scaling of two temperature correlators. We close 
by concluding that C has a value compatible with our 
generalized universality picture but not with the naive 
application of Eq. ^ (be. assuming T = values for the 
domain- wall and chaos exponents). 

Model, observables and methods - We analyze the 
2D Edwards- Anderson spin glass with Ising spins (o~i — 
±1). Its Hamiltonian has the form 

Hj{°) = ~Y1 JijO-iO-j , (2) 

<ij> 

where the sum runs over all pairs of nearest neighbor 
sites on a square 2D lattice of linear size L with periodic 
boundary conditions. The couplings are independent 
quenched random variables taking the two values ±1, 
each with probability 1/2. 

To analyze the possible chaotic behavior of the system, 
we look at correlation functions of observables computed 
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at two different temperatures. For a generic observable 
O we use the correlation function 



Co(L, T, AT) ee {<D[T) ~ °( T ))(°( r ') - °( T ')) 



\j (O(T) - 0(T)) 2 \f (O(T') - 0{T')f 



(3) 

where T' = T + AT, L is the linear size of the lattice 
and the bar stands for the quenched average over the 
distribution of the couplings Jy . In the rest of this note 
we will focus on two main cases: in the first one O is the 
total free energy F of the system, while in the second 
case it is the domain wall free energy, Ff,w- 

We consider a L x L lattice with given quenched 
couplings Jij, and we determine F(T) and Fdw(T) 
via the exact computation of the partition function 
Z ^llHjEl- We can then average over different samples 
to obtain (up to some statistical errors) the correlation 
functions C (L,T,AT). 

The total free energy - We start by discussing the 
computation of the free energy, where we identify O with 
the total free energy F of the system, and we compute 
Cf{L,T, AT). In a system where the total free energy 
has a chaotic behavior we would expect that in the large 
L limit Cf should drop very fast [2fJ, as a function of 
L, even for infinitesimal values of AT. In the following 
we analyze Cf and we do not find a chaotic behavior. 
We use these results to show that even in the T — ► 
limit the total free energy does not behave chaotically in 
temperature. 

We have averaged over a number of samples for differ- 
ent values of the system size V ee L x L. Specifically, we 
have collected 99809 samples at L = 10, 31571 samples at 
L = 20, 3117 samples at L — 30, 818 samples at L = 40 
and 336 samples at L = 50. Our error bars have been 
computed by using the jack-knife method (for details see 
21]). 

In our approach we obtain the full density of states, 
that allows us to compute expectation values at every 
temperature (or couple of temperatures) we are inter- 
ested in. Because of that we select our temperature 
ranges (both for T and for AT) using the physical be- 
havior of the system: we want to be for example as much 
as possible in a scaling regime, i.e. at low T, but not at 
very low T (i.e. where, at given L, the system al- 
ready shows an unphysical behavior dictated by the pres- 
ence of a gap). Because of these facts we will present in 
our study data that span a range of temperatures going 
from T = T min ~ 0.2 to T = T max ~ 0.4: these val- 
ues are chosen such that, for the considered linear sizes, 
we are in the scaling region where the relevant physical 
phenomena occur. 

In figure we plot In Cp (L, T, AT) as a function of 
AT, for T = 0.25; (that we know from [Til Hg| to be 
in the scaling regime for L ~ 50); we use in the plot 
a AT step of 0.01 (reasonably small with respect to T). 
Somehow already this first, simple plot hints that the free 
energy of the system does not unveil a chaotic behavior. 
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FIG. 1: log{CV(T T, AT)} versus AT, for T = 0.25. 



In fact one sees that the correlation function rises when 
L is increased. This trend is even more evident for larger 
values of AT. Such a behavior is the opposite of what 
we would expect in a chaotic scenario, where for large L 
values no correlation survives. 

The issue of a possible presence of chaos in the free 
energy in the limit T — > (i.e. when T approaches the 
T = critical point) has to be dealt more carefully. For 
doing that we need a quantitative analysis of our data. 
We start by taking the infinite volume limit of our cor- 
relation functions. We present here the analysis of the 
correlation function with T = 0.25, but the other T val- 
ues follow the same pattern. We have analyzed a number 
of temperatures in the scaling region that we are able to 
access (say for T going from 0.2 up to 0.4 [lfil llfll |). For 
each value of AT we fit C F (L, T = 0.25, AT) to the func- 
tional form: 



C F {L,T = 0.25, AT) 



C'f(T = 0.25, AT) 
A(T = 0.25, AT) 
L ' 



(4) 



determining in this way the large L limit C f (T, AT) . 
As we noticed already, the correlation function increases 
with L: the fit to the form of Eq. 01 works very well, and 
Cf reaches a strictly positive asymptotic value in the 
thermodynamic limit. 

From there, we are able to study C F {T, AT) as a func- 
tion of AT (and, again, we do this for different values of 
T). A very good fit is obtained by assuming the func- 
tional dependence: 



C F (T, AT) = exp (-a(T)AT 2 ) 



(5) 



In all cases this fit works very well; for example at T = 
0.25 we estimate a ~ 0.114. We find that in the limit 
T — > the function a(T) smoothly goes to a constant 
value (as opposed to its possible divergence, that would 
indicate chaos in the free energy): in the limit T — * 
there is no chaos in the total free energy of the system. 
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Domain walls and their fractal dimension — Our sec- 
ond observable is based on domain walls. We introduce 
a domain wall into the system by applying to a given 
realization of the random quenched couplings first peri- 
odic boundary conditions (pbc) and then anti-periodic 
boundary conditions (apbc). For each sample, at any 
temperature, we define the domain-wall free energy as 
the difference of the free energies of these two systems; 
the same can be done for the energy and entropy. So for 
each realization of the couplings we compute the parti- 
tion function twice: once with pbc and once with apbc 
(equivalently, we compute the partition function of two 
different systems with pbc, where in the second system 
we have inverted the sign of one line of couplings of the 
first). 

Before discussing our results for this interesting corre- 
lation function we discuss a byproduct of this approach, 
that will be crucial in our final discussion of the physics of 
two-dimensional spin glasses. A nice feature of the exact 
partition function approach is that it gives us in particu- 
lar the entropy of the domain wall even at zero tempera- 
ture: this is a quantity that vanishes in the models with 
continuous Jij coupling values, but that is non-zero in the 
model with binary quenched random couplings where the 
ground state degeneracy is high. 

The zero temperature domain-wall entropy ASdw has 
a disorder mean of zero (since inverting a number of cou- 
plings maps us to another sample in the disorder ensem- 
ble that has the same probability but has the opposite 
value of ASdw)- We thus focus on the mean square 
fluctuations; these have a characteristic size growing as 

\ASdw\ ~ irt (6) 

when L — > oo. Previous work by Kardar and Saul |l3T ] 
(also based on partition function computations, but with 
smaller statistics and for smaller L values than us) gave 
the estimation d s « 1.0. We have used our data to obtain 
a more precise measurement of d s . In Fig. [2 we show the 
scaling with L of the absolute value of ASdw (to exhibit 
the very accurate power scaling) . We also show our best 
fit which gives d s = 1.03 ±0.02 (where we only quote the 
statistical error). This exponent d s is usually interpreted 
as the fractal dimension of the domain wall created in the 
system; when d s > 1, the domain wall is rough. However, 
in the ±J model there are many ground states and so it 
is not appropriate to think of d s as associated with one 
interface; in fact it is not a priori necessary that d s > 
1. We note again that for Gaussian quenched random 
couplings one finds d s = 1.27 [ill H^: we thus see that 
having many domain walls in our degenerate system leads 
to smaller entropy fluctuations. 

The chaos exponent £ - Let us now discuss temper- 
ature chaos in the domain wall free energy Fdw, defined 
as the difference of the free energy in the system when 
using periodic boundary conditions versus anti-periodic 
boundary conditions. Given this quantity we use the 
definition © and we compute Cf dw - Notice that the 
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FIG. 2: \ASdw\ versus L. 
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FIG. 3: log{CV DK - {L, T, AT)} versus AT, for T = 0.25. 



relation © is simplified from the fact that Fdw = 0. 

Here we have used slightly smaller sizes than in the 
other case. For the linear sizes L = 10, 12, 20, 24, 
30 we use different numbers of samples: 99808, 94351, 
31570, 4098, 3116 respectively. In figure |3| we plot 
lnC FDW (L, T, AT) as a function of AT for T = 0.25. 

Fig. |3 is very different from Fig.^ the linear correla- 
tion coefficient of the domain wall free energies at two dif- 
ferent temperatures (with AT being fixed) decreases as L 
increases. As opposed to the case of the total free energy 
of the system, this trend is compatible with an emergent 
chaotic behavior. The same pattern arises when looking 
at the domain wall energy (Edw = E p i, c — E ap b c ) and at 
the domain wall entropy, (Sdw = Spbc — S ap bc) and ana- 
lyzing the related correlation functions Ce dw (L, T, AT) 
and Cs DW (T, T, AT). 

The fact that the total free energy does not show any 
sign of chaos while the domain wall free energy hints for 
a possible chaotic behavior is, as we have already dis- 
cussed, very natural. The small value of the exponent 9 
implies that in the domain wall free energy large cancella- 
tions play an important role |2| |3j . This suggests in turn 
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that the domain wall energy and the domain wall entropy 
(times T) exhibit a very similar temperature dependence 
and greatly contribute to this cancellation. This behavior 
has as a natural consequence the presence of temperature 
chaos: as soon as the parameters of the system get mod- 
ified, even very slightly, cancellations will act in a differ- 
ent way and select completely different domain walls (the 
same arguments suggests the presence of chaos under a 
small change in the quenched disorder). 

In the scaling theories, the domain wall free energy 
scales with the lattice size as L 6 (where 9 is the stiffness 
exponent) and the domain wall entropy scales as L d ^ 2 
(where d s is the fractal surface of domain wall) Q . Con- 
sider first the system at temperature T. Here the domain 
wall free energy scales with L as: 



Fdw(T) = E DW — TS 



DW 



Y(T)L 6 



(7) 



where Y(T) is the generalized stiffness coefficient. Now 
if the system is at temperature (T + AT), the domain 
wall free energy scales as: 



DW 



(T + AT) = E DW - (T + AT)S DW 
« Y{T)L 6 - ATL d °/ 2 . 



(8) 



Since the droplet theory predicts that d s /2 > 8, the quan- 
tity (JSJ for large L can change the sign between T and 
T + AT. This implies that the equilibrium state changes 
when one reaches a length £ such that 



that is 



AT £ ds/2 



1 \ d s /2-e _ / 1 



AT 



f 

VAT 



(9) 



(10) 



It is important to note that these scaling laws are valid 
for small T > 0, and so both and d s should be construed 
as being obtained at T > 0; this subtlety will be essential 
for interpreting our data. 

When AT is small and all relevant length scales diverge 
one expects that |5j 



C Fdw (L,T,AT) 



L 



£{T, AT) 



2C 



(11) 



for L smaller than £(T,AT), i.e. for small enough AT. 
In other words in this limit I — Cf dw (L,T, AT) scales 
as a power of L. In figure 01 we plot (I — Cf dw ) vs. L 
on log-log scale for small values of AT (AT = 0.01 , 0.03, 
0.05, 0.07), when T = 0.25 (similar results are obtained 
for other values of T). The scaling in L allows us to 
obtain £ and £(T, AT) from a linear fit to ln(l — Cp nw ); 
we obtain ( values between 0.92 and 0.99, depending on 
the values of T and AT. 

We have also extracted from our fits the quantity 
£(T,AT). Fig. [5] shows that this length scale diverges 
as predicted in Eq. The different best fit values for 



0.01 



0.001 



0.0001 




FIG 

T = 



4: (1 — Cf dw ) vs L in double logarithmic scale. Here 
0.25. 




1/AT 



FIG. 5: £(T, AT) versus -^p in double logarithmic scale. 



£ can be summarized by quoting ( = 0.95 ± 0.05; we do 
not detect any systematic dependence on T. 

In Fig. we show an interesting data collapse. We 
plot 1 - C Fdw as a function of L/£(T,AT), for T = 
0.25 and T = 0.3 in double logarithmic scale. The fact 
that we get a single curve shows that the assumption 
that Cf dw (T, T, AT) is only a function of L/£(T, AT) 
is reasonable. The fact that the curve is straight gives 
further support to Eq. 1111 

Discussion and conclusions - This note spells a 
number of results. Let us start by summarizing them. 
First, when considering 2D Ising spin glasses with binary 
couplings, the free energy F does not show any tempera- 
ture chaotic behavior, even in the T — * limit. This is in 
contrast |20| with what happens in the directed polymer 
in a random medium, but in line with the prediction [(| 
of the Migdal-Kadanoff approximation. A simple inter- 
pretation is that the sample to sample non-chaotic fluc- 
tuations of F are O(L), far larger than the chaotic O(L ) 
fluctuations arising from droplets. 

As opposed to the total free energy, the domain wall 
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free energy does behave chaotically in T. The question of 
interest is whether the ± J model has chaotic behavior in 
the same class as the continuous distribution models. We 
first measured the fractal dimension of domain walls at 
T = exactly and found dj =0 ~ 1.0. Then we measured 
directly the temperature chaos exponent £ in the limit of 
small T and found ( = 0.95 ±0.05. (Note that these two 
measurements are conceptually very different.) 

The point of view introduced by ^5| and |l6j is the 
following. If T > 0, all 2D Ising spin glasses, indepen- 
dently of the distributions of the quenched couplings, 
enjoy a strong, generalized universality a binary 

distribution of the couplings gives rise to the same criti- 
cal behavior as the Gaussian distribution. Only exactly 
at T — a very peculiar set of coupling distributions 
(including binary couplings) generates an anomalous be- 
havior 15] : for cases where, by composing an arbitrary 
number of couplings, an energy gap survives (for example 
a distribution where J = ±2, ±3 is in this class, while a 
distribution where J = ±\/2, ±V3 is not), then at T = 
one gets 9 = 0. Our findin gs h ere for T — ► complete 
the scenario designed in [lfitllfij and are completely con- 
sistent with it. If we use indeed the relation implied by 
the droplet theory £ = d s /2 — 9 and plug in both the 
value 9 ~ —.285 of the model with Gaussian quenched 
couplings and the value d s ~ 1.27 obtained by [l4| and 
by [1 1| we would expect £ ~ 0.92, in very good agreement 
with our estimate £ T ~*° = 0.95 ± 0.05. We believe that 
this is a very important point of support of the picture 
proposed and advocated in ^5|, [lif: there is one single 
critical theory that emerges, when T — > in 2D spin 
glasses. The exponents of this theory are, at this point, 
well determined. 

Our measurements also confirm that T — is for the 
,/ = ±1 distribution, a singular point. In that model 



at T = 0, we found d%=° ~ i.o. This, together with 
the value 9 = 0, implies that, exactly at T = 0, if the 
relation £ = d s /2 — 9 were to be valid, we would have 
£ T =° = 1/2. The difference of the physics at T = and 
what one observes when T — + (as defined from the order 
used for taking the T — ► and the L — > oo limits) is at 
this point self-evident. 
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